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SUMMARY 


A method is presented hy which the pressure distribution on a 
slender, arbitrary "body of revolution moving at supersonic speeds at 
zero angle of attack may be calculated with an appreciable saving of 
time as compared with that required by other methods. The method was 
developed from the linearized potential— flow relations given by Jones 
and Margolis in NACA TN No. 108l and is essentially a graphical computing 
technique based on the assumption that a body may be represented satis- 
factorily by a relatively small, number of point sources and sinks. An 
example is presented to illustrate the computational procedure. 

The method gives results for parabolic bodies in excellent agreement 
with results obtained by the method of Jones and Margolis. It was found 
that a parabolic body of fineness ratio 6 could be represented satis- 
factorily by as few as 10 sources and, sinks. 

Pressure distributions computed for the German missile AWlP 
at Mach numbers of I .87 and. 1 . 56 are compared with wind-tunnel measure- 
ments and with results of calculations by the characteristics methods 
of Sauer and Tollmein to provide some information on the accuracy and 
reliability of the method presented herein when applied to arbitrary 
shapes. Although the body of the A^J-VIP is somewhat more blunt than 
would be desirable for treatment by linearized methods, the method 
presented reproduces the salient features of the measured pressure 
distribution with comparatively small, discrepancies. 


INTRODUCTION 


The problem of computing the pressure distribution on a body of 
revolution moving at supersonic speeds has been treated by several 
authors (references 1 to 6) with varying degrees of appro xi ma t ion . 
References 1 to 5 present methods based on the isentropic linearized 
potential— flow equation first applied to the problem by Von Harman 
in reference 1. A more rigorous method based on the method of character- 
istics is also presented in reference 5 and is generalized by Ferri 
in reference 6 to include nonisentropic rotational flows. The usefulness 
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of these methods for the design of supersonic aircraft or missiles has 
been severely limited, however, by the excessive labor involved in 
the computations . 

Jones and -Margolis (reference 4) • and Lighthill (reference 3) have 
shown that, for certain mathematical body shapes, the pressure distribution 
may be obtained in the form of relatively simple equations. In general, 
however, the designer cannot utilize these simple mathematical shapes 
exclusively because of other requirements of the design; for example, 
providing vision for the pilot or seeker, tankage space equally distributed 
about the center of gravity, and shape changes to allow simpler or 
stronger construction. In addition, the possibility exists that body 
shapes other than those specified by simple mathematical relations might 
have the lowest drag, particularly when these • shapes are combined with 
wings to form a practical airplane or missile configuration. The need 
of a rapid method for computing the pressure distribution on arbitrary 
bodies is therefore evident. 

An approximate method for computing the pressure distribution on 
slender, arbitrary bodies of revolution moving at supersonic speeds and 
zero angle of attack is presented herein. Use of a graphical computing 
technique allows the pressure distribution to be obtained at an appreciable 
saving in time as compared with other methods. The method (which is 
referred to hereinafter as the "point-source method") was developed from 
the results of Jones and Margolis (reference 4) and Is based on a 
simplifying assumption concerning the source distribution used to represent 
the body. The validity of this assumption in the case of parabolic bodies 
is e xa m ine d by comparison of results obtained by the point-source method 
with results obtained by the method of Jones and Margolis. 

Although results obtained by the point-source method are subject to 
the uncerta int ies inherent in the linearized potential— flow equations 
and to additional uncertainties due to the' graphical approximations used 
in the computations, the method is believed to provide a quick and useful 
means for comparing the supersonic characteristics of slender, arbitrary 
bodies of revolution. Some information concerning the reliability of 
the results was obtained by comparison of pressure distributions computed 
by the point-source method with calculations by the characteristics method 
and with wind— t unn el measurements presented in references 7 and 8 for the 
German missile A^YIP. 

The computational procedure for the point— source method is 
illustrated by an example . 


SYMBOIS 


P pressure coefficient 

x,r 


cylindrical coordinates 
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| abscissa of source or sink (measured along x-axis) 

f 1 (| ) rate of change of source strength (with x) 

R radius of body 

F fineness ratio of body 

9 half, angle of cone 

Q pressure function (f t (|)Ax) 

Ax: interval length 

n positive integer 

k constant 

M Mach number 


0 = \/M 2 - 1 


Subscripts indicate that the value of the quantity being considered 
is taken at that particular point. Subscript o refers to the nose of 
the body. 


DEVELOPMENT OF METHOD 


In reference 4 Jones and Margolis show that a slender body of 
revolution moving at supersonic speed with zero angle of attack may be 
approximated by an infinite number of sources or sinks located on the 
axis of the body and, that the pressure coefficient P at any point x 
on the body surface is given by 


ix-PRx 


P x = 2 


f*U) d£ 


1*0 


)/(x- i) 2 - (m x ) 2 


(1) 


where f * ( |) is the first derivative of the function 


f(i) = 


b m 


(2) 




Explicit relations may be .easily obtained from, equations (l) and (2) for 
the pressure distribution -on cones and on parabolic bodies. For arbitrary 
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"bodies., the pressure distribution can also be computed directly from 
equations (l) and (2) (provided that f(|) is a continuous function 
of x) by either analytical or graphical procedures; however* the 
graphical approach involves extensive computation and a separate graphical 
integration for each point considered* and the analytical approach 
involves the approximation of the body shape by a suitable expression 
and (except in special cases) difficult and tedious integration. 

Equation (l) may be written as 



2Qi 


x - £ ± ) 2 - (PR*) 


2 


(3) 


where each of the n point sources located on the axis of the body 
ahead of (x — PR X ) represents a section of the body of length Ax. 
Comparison of equation (3) with equation (l) shows that 


Q = f*U)Ax 




It is evident that* if the arbitrary body may be adequately represented 
by a relatively small number of point sources* ^ the pressure distribution 
may be calculated rapidly by use of equation (3). The computation is 
greatly facilitated by the direct proportionality between the pressure 
coefficient and the value of Q for a single source. (See equation (3) 
for n = 1.) This direct proportionality a l lows the pressure field for 
a single source to be computed once as. a function of Pr and used for 
any source or sink: at any Mach number within the lim i tations of the 
basic equations. This universal pressure field may be conveniently 
plotted as a "grid" of lines of constant pressure coefficient on 
the x*Pr plane. Values for this grid computed from equation (3) 

(for n = 1* Q = 0.01* | = 0; that is* a single source having a convenient 
value of Q and located at the origin of the grid) are presented 
as table I. The grid is plotted in figure 1. The pressure distribution 
over the body due to a single source representing any section of the 
body of length Ax may be readily obtained graphically by superimposing 
the origin of the grid at the source location on a drawing of the body 
shape (plotted as PR against x)* and reading the pressure shown by 
the grid at a number of points on the body surface* and then multiplying 
these pressures by the ratio of the Q for the source being considered 
to the value of Q for which the grid was prepared. The process is 
repeated for each source necessary to represent the body* and the pressure 
coefficient at any point is obtained by adding the pressures at the 
point due to each source. 


If the body is to be satisfactorily represented by a s m a l l number 
of discrete sources* the location of a source within its section of the 
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"body will "be Import ant . It was found (see appendix) that the point source 
should he located on the axis Ax/4 ahead of the rear end of the section 
of the body of length Ax which it represented and the Q for the 
source should be computed from the radius and slope of the body surface 
at the source location. With this source location and strength, the 
pressures on the body surface at and behind, the intersection of the Mach 
line from the axis at the end of the interval with the body surface 
closely approximated the pressures obtained from a cont inuo us distribution 
of sources along the axis of the section. As the pressures on the body 
due to a number of sources are summed only at the intersection of the 
Mach line from the end of each interval with the body surface (see 
appendix), the infinite pressures which equation ( 3 ) gives for x = | + pr 
(along the Mach line from each source) are avoided. The validity of 
results obtained by the point— source method and the computational 
procedures used are demonstrated by an example. 


EXAMPLE 


Computation of a pressure distribution by the point— source method 
is accomplished in the following steps: 

(1) Divide selected body shape into sections and determine source 
locations . 

(2) Graphically compute pressures on body due to sources of Q = 0.01 
by successive superposition of grid on body shape plotted as PE against 

( 3 ) Prepare tabular form for .recording values. 

(4) Determine value of Q for each source. 

( 5 ) Compute pressure on body due to each source. 

( 6 ) Sum pressures due to each source to obtain desired pressure 
distribution. 

A parabolic body of fineness ratio 12 was chosen for this example 
as representative of bodies suitable for treatment by the linearized 
theory and because it could be easily treated by the method of Jones 
and Margolis (reference 4) for comparison. The pressure distribution 
was computed at a Mach number of 1.4l. 

Body.— The body shape is specified by 

E = i(l - 4x 2 ) ( 5 ) 

Equation ( 5 ) describes a parabolic body of unit length symmetrical about 
its maximum diameter which is located at x = 0. The body was arbitrarily 
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divided into 20 sections of equal length (Ax = 0.05), and, an described 
in the preceding section, a source location was chosen Ax/4 ahead of 
the end of each section. 

Graphical computations .— The successive superposition of the grid 
origin at the chosen values of £ is accomplished by plotting the grid 
fl-nrl body shape to the same scale on separate sheets of semitransparent ' 
coord inat e paper and matching the two sheets on a glass surface lighted 
from behind. A convenient size for the grid for most purposes has been 
found to be 10 by 15 inches with the Pr— scale 10 times the x— scale. 

The grid (fig. l) should be plotted carefully from the data of table I 
and can, of course, be used for a] 1 bodies and Mach numbers. The body 
shape (plotted in the form 0E against x) is shown in figure 2 with 
the grid origin superimposed on the source £ = — O.2625 which represents 
the section of the body between x = -0.30 and x = -O.25. 

Table.— A convenient form for tabulating the values obtained in the 
computation is presented as table H. The abscissas of each of the 
20 sources are arranged vertically in two sets on the left of the table 
and the values of P used for the grid are used for column headings. 

The values of x at which the lines of constant P of the grid cross 
the body surface are recorded in the appropriate column in the upper part 
of the table. For the case shown in figure 2 (£ = -0.2625) the values 
of x in table H are underlined. 

The last column in the upper part of table H, headed x^, gives the 
values of x at which the pressures due to each source are to be summed. 
These values are obtained by superimposing the grid origin at the end 
of the interval end reading the value of x on the body where the 
line of P = 00 (Mach line) crosses the body surface. For the interval 
just mentioned (| = -O.2625), the grid origin is set at x = -0.25 a 11 ^- 

the value found for x^. is -0.216. 

Determination of Q for each source .— In this example the values 
of Q for each of the 20 sources may be found easily by use of 
equations (2), (4), and (5). In the more general case of the arbitrary 
body with which this paper is concerned, however, the value of Q for 
each source must be dete rmin ed from the ordinates of the body by use 
of the differential forms of equations (2) and (4) as shown in figure 3* 
Any convenient incr eme nt may be used to determine the variation 
of B(dR/dx) with x from the variation of E with x (lower part 
of fig. (3))* However, as the value of Q must be deter mine d at the 
source location, A[R(cLR/dx)] is taken from the midpoint to the end of 
the interval, and multiplied by 2 (upper part of fig. 3) • The ratios 
of the value of Q dete rmine d for each source to the value of Q used 
to prepare the grid (O.Ol) are tabulated in the last column of the lower 
part of table II opposite the value of 5 T° r each source. 

Computation of pressures on body due to each source .— The data in the 
upper part of table H is the pressure distribution over the body if all 
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the sources have Q = 0.01. As the pressure at any point due to one 
source is directly proportional, to the value of Q, the actual pressure 
distribution on the body due to the source is determined by multiplying 


the value of P heading each column by the ratio — for the 

Qgrid 

particular source being considered. The resulting P is entered in the 
lower part of table U . For example, for the source i = — O.2625, the 


x = -0.210 


ratio — is -0.0118 and the pressure coefficient at 

Qgrid 

(upper part of table, line | = — 0.2625, column P = O.5) is 
P = -0.0118(0.5) = -0.0059 (lower part of table, line £ = -O.2625, 
column P = 0.5)- For the same source the pressure coefficient 
at x = — 0.132 (upper part, line | = —0.2625, column P = 0i.l6) 
is P = -0.01l8(0.l6), = -0.0019 (lower part, line | = -O.2625, 
column P = 0.l6). These values are bracketed in table II. 


Summation of pressures due to each source .— The pressure coefficients 
for each source tabulated in the lower part of table II are plotted against 
the corresponding values of x at which they occur (upper part of table II) 
in figure 4. The logarithmic ordinate scale is used so that the points 
are separated at low values of the pressure coefficient. Curves are 
faired through the points for each source and all of the curves are added 
at the summation points previously determined (last column, upper part 
of table H). These sums and the values of x^ at which they occur are 
the desired pressure distribution on the body. The pressure coefficient 
at the nose of the body is computed directly from equation (Al) 

(see appendix) under the assumption that the extreme nose of the body 
approximates a cone. 


Computation for different Mach numbers .— In order to compute the 
pressure distribution on the same body at a different Mach number, the 
body shape (in the form PR against x) must be replotted and the data 
in the upper half of table II redetermined. The data in the lower half 
of table II are not changed and the pressure distribution is computed 
from the new table exactly as in the preceding example. 


Time required for computation. — Once the grid has been prepared, 
computation of a 20-point pressure distribution for an arbitrary body 
by this method required about 8 manhours and, for the Same body, results 
may be obtained for a different Mach number in about -4- additional man- 
hours. For parabolic bodies only the same number of points can be 
computed by the method of JoneB and Margolis (reference 4) in somewhat 
less than half this time. 
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EVALUATION CSF METHOD 


Comparison with results by method of Jones and Margolle .— The 
pressure distribution determ in ed. In the example (fig. 4) Is reproduced 
in figure 5 together with the pressure distribution computed for the 
same body by the method of Jones and Margolis (reference 4) . The 
expression for this pressure distribution may be obtained from 
equations (l), (2), and ( 5 ) and is (see reference 4) 



-1 


+ 6 [(x - xj - lj-x] \/(x - Xq) 2 - (pRx) 2 f 


(6) 


E x a m i n ation of figure 5 shows that the pressure distribution determined 
by the present method agrees with that of Jones and Margolis except for 
a small and apparently random scatter which approximately corresponds 
to the uncertainty of the graphical and numerical computations. 

In order to investigate further the uncertainties introduced by the 
approximations involved in the point— source method, the pressure distri- 
bution was calculated for a parabolic body of fineness ratio 6. 

Inasmuch as large pressure coefficients and gradients are encountered 
for the body of fineness ratio 6 at M = 1.4l, discrepancies should 
be more easily discernable than in the case of the body of fineness 
ratio 12. The results of this calculation for 10 and 20 intervals 
(Ax =0.10 and O.O 5 ) are compared in figure 6 with results calculated 
from equation (6) . This fineness ratio was chosen because it approxi- 
mates the lower limit of fineness ratio for which the linearized theory 
gives a completely supersonic flow field about the body at a Mach 
number of l.hl. It is interesting to note, however, that recent results 
obtained by the NACA wing— flow method showed that pressure distributions 
of the type predicted by Jones arid Margolis occurred for a body of 
fineness ratio 6 at a much lower supersonic Mach number even though a 
large subsonic region was present at the nose. 

Figure 6 shows that once more the results of the present method 
agree closely with those of Jones and Margolis. The scatter of the data 
in figure 6 (body of fineness ratio 6) appears to be somewhat greater 
than the scatter of the data in figure 5 (body of fineness ratio 12) . 
Such a result might be expected because of the steeper gradients 
encountered on the body of fineness ratio 6. Little difference is 
evident in figure 6 between the results for the two different intervals. 
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As additional computations for the intervals of 0.02 (50 sections) near 
the nose of the body showed no significant difference from the results 
presented, it is concluded that the results obtained are not appreciably 
affected by the interval size (within the range investigated) and that 
a parabolic body may be satisfactorily represented by as few as 10 sources. 
The method does not reguire that a constant interval length be used and 
for arbitrary bodies it is recommended that shorter intervals be used in 
regions where the body cross section is changing rapidly. 

✓ 

The small discrepancies between the results computed by the point- 
source method and the method of reference 4 are negligible in comparison 
with the uncertainties introduced into both methods by the assumptions 
involved in the linear theary from which equations (l) and (2) are taken. 

The linear— theory results have been shown by several authors (for 
example, reference 2) to be in error by as much as 20 percent at some 
Mach numbers in the case of a 10° cone. 

Comparison of results by point— source method with results by 
characteristics methods and with wind-tunnel measurements .— Some 
information on the usefulness and reliability of the point-source method 
may be obtained by comparison of results with results of calculations 
by the characteristics method and with wind-tunnel measurements for 
the German missile AWlP presented in references 7 and 8. Pressure 
distributions measured at Mach numbers of I.56 and I.87 are presented 
in reference 7 and these measurements are compared with calculations by 
the graphical characteristics method of Tol Imein and the computational 
characteristics method of Sauer in reference 8. The body shape of 
the AWlP is not particularly suited for calculation by linearized 
methods because of the bluntness of the nose (ogival nose having a nose- 
length to maximum diameter ratio of 3 - 9 ) nose half angle 15 • 9° - however, 
the AkVlP is the only body for which relatively complete experimental 
data are available, and the comparison should indicate the accuraoy with 
which the method predicts the pressure distribution on such relatively 
blunt bodies. 

The measured pressure-distribution points and fairing (solid line) 
reproduced from reference 7 are shown in figure 7 together with the . , 
results obtained by the method presented herein (dash line). Examination 
of this figure shows that the point— sourcej.method satisfactorily reproduces 
the salient features of the measured pressure distribution but shifts the 
sharp breaks which occur at the beginning and end of the cylindrical 
section toward the rear of the body. The magnitude of the maximum positive 
and negative pressures are also somewhat exaggerated. 

The fairing of the experimental data and the pressure distribution 
calculated by the 'point— source method shown in figure 7 a* 1 © reproduced 
in figure 8 where they are compared with calculations by the characteristics 
methods of Sauer and Tollmein (reference 8). Figure 8 shows that Sauer’s 
calculations agree closely with the experimental results and Tollmein* s 
agree less closely, particularly near the middle of the nose and on the 
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cylindrical section of the body. The calculation by the point-source 
method agrees closely with that of Tollmein near the nose of the body, the 
most ohwious discrepancy being the exaggeration «nd rearward shift of 
the minimum pressure peak near the beginning of the cylindrical section 
previously mentioned. 

A comparison similar to that shown in figure 8 is presented in 
figure 9 for a Mach number of 1.56. In general, this comparison shows 
the same results as does figure 8; however, the overestimation of the 
pressures on the nose by the characteristics methods, which was barely 
visible in figure 8, is increased appreciably at the lower Mach number. 

\ 

It will be noted that none of the calculated pressure distributions 
agree with the measurements on the rear ogive of the missile. Beference 7 
states' that the pressures in this region are greatly affected by the 
four large stabilizing tail fins which were in place when the measurements 
were made. Also, these pressures could be subject to large boundary- 
layer effects. 

The sharp breaks in the pressure distributions for the A4Y1P are 
characteristic of the class of bodies having cylindrical center sections 
faired into a curved nose and/or tail section so that the second 
derivative of the radius with respect to x is discontinuous. It is 
apparent, however, that the minimum pressure peak computed from 
equations (l) and (2) always occurs at a distance PE behind the 
beginning of the cylindrical section rather than at the beginning of the 
cylindrical section where the experimental results and rigorous theoretical 
treatment show it must be located. The occurrence of the rearward shift 
for one class of bodies suggests that such a shift may be present in the 
pressure— distribution results computed from equations (l) and (2) for 
other types of bodies. Extensive experimental pressure— distribution 
measurements are required to verify the existence of the shift for 
arbitrary bodies, however. 

In a sm uch as the approximations involved in the linearized theory 
become more representative as the body becomes more slender and because, 
in general, slender bodies are better suited for aerodynamic purposes, 
the method presented herein is believed to provide a useful means for 
rapid estimation of the pressure distribution on slender, arbitrary 
bodies of revolution suitable for use as supersonic aircraft or missiles. 
Insufficient experimental evidence is available to establish definitely 
the limits of applicability of the method; however, it appears doubtful 
that results of usable accuracy could be obtained for bodies more blunt 
than the A4Y1P used as an example herein. 

Illustration of change in pressure distribution resulting from a 
moderate change in body shape .— In order to illustrate the magnitude of 
the changes in pressure distribution which result from moderate change 
in body shape the pressure distribution calculated for the test body of 
fineness ratio 12 used in the Langley Eree-Eall Test Program is shown 



MCA TN No . 1768 


11 


in figure 10 compared with that for a parabolic body of the same fineness 
ratio. The body shapes are compared in the lower part of the figure. 

The free-fall test body was derived from an MCA low-drag airfoil and 
had a cusped tail faired into a tail boom and a relatively blunt nose. 

Some indication of the over-all accuracy of the pressure distribution 
calculated for the free-fall test body may be obtained by comparison 
of the drag coefficient obtained by integrating the pressure distribution 
with the experimentally measured value. The computed value (0.080) is 
in good agreement with experimental results if a skinr-friction coeffi- 
cient of 0.003 based on surface area is assumed. Examination of 
figure 10 shows that the free-fall test body is more critical than the 
parabolic body with respect to flow separation because of the steep 
adverse gradient encountered near the rear of the body. 


CONCLUDING- REMAKES ' 


A method is presented by which the pressure distribution on an 
arbitrary body of revolution moving at supersonic speeds may be calculated 
rapidly. The method was developed from the linearized potential— flow 
relations reported by Jones and Margolis in MCA TN No. 108l and by 
other authors. 

The method is essentially a graphical computing technique based on 
the assumption that a body may be represented satisfactorily by a 
relatively nmnl 1 number of point sources and sinks. RreBsure distri- 
butions computed for parabolic bodies by this method were in excellent 
agreement with the results of Jones and Margolis . As few as 10 sources 
and sinks satisfactorily represented a parabolic body of fineness 
ratio 6 . 

Comparison of pressure distributions calculated for the German 
missile AtVlP with wind— tunnel .measurements showed that the point— source 
method reproduced the salient features of the measured pressure 
distribution, but that the minimum pressure peak was exaggerated and 
shifted to the rear. 

Pressure distributions for the AkYU? at Mach numbers of I .87 and I .56 
computed by the point— source method were somewhat less satisfactory 
than those computed by the characteristics methods of To llme in or Sauer. 
However, as the time required to compute a pressure distribution by 
the point— source method is appreciably less t han is required by either 
of the two characteristics methods and as the salient features of the 
pressure distribution are satisfactorily reproduced, the point-source 
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method is believed to provide a useful means for rapid estimation of 
the pressure distribution on slender, arbitrary bodies of revolution 
suitable for supersonic aircraft or missiles. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Ya., September 23, 19^ 
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APPENDIX 

LOCATION OF A SINGLE SOURCE EQUIVALENT TO A CONTINUOUS 
DISTRIBUTION OF SOURCES WITHIN A FINITE INTERVAL 


In the linearized theory of supersonic flows (references 1 to 5 ) 
bodies of revolution are assumed to be represented by a continuous 
distribution of sources and sinks along the axis of the body and the 
pressures resulting from these sources and sinks are assumed to be 
propagated only within the Mach cones. Thus, if an interval of 
length Ax extending from the Vertex of a semi— Infini te cone is 
considered, it is apparent that only at the intersection point (the - 
intersection of the Mach cone from the axis at the end of the interval 
with the cone surface) is the pressure influenced by al 1 of and only by 
that part of the continuous source distribution within the interval - As 
the rate of change of source strength f* (g) is constant for a cone, 
the location of a single source which causes the same pressure at the 
intersection point as does the continuous distribution of sources within 
the interval may be determined in the following maimer. 

pressure coefficient on the surface of a cone ; of semivertex 
due to a continuous distribution of sources along the axis may 
from equations (l) and. (2) to be constant and equal to 

P = 2tan% cosh -1 ^—-^ (Al) 

pressure coefficient on the surface of the cone at any point x 
single source located at i is (see equation ( 3 ) for n = l) 

B ' 2Q 

** ~ / ■■ = - (A2). 

/U - l ) 2 - (PR*) 2 

As for the cone (with vertex at the origin) 


The 
angle Q 
be found 


The 
due to a 


R = x tan 6 


then from equations (2) and (^) 


Q = Ax tan % 
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Equating the pressure coefficients given "by equations (Al) and (A2) 
at the intersection point and expressing the location of the source ahead 
of the rear end of the interval as a fraction of the interval length 

(& = Ax - t, r = - C ? . 4 J ^ gives 



Computations from equation (A3) for values of r from 1.01 to 10.0 
show that the ratio 5 /Ax is approximately equal to O.25. The maximum 
variation in this range , -which corresponds to cdne angles of 5° to 45° 
at M = 1.4l (or Mach numbers -Trom 1.15 to 5*75 for a 10° cone), is 
less than ±1 percent of Ax. 

An expression for the pressures on the surface of the cone behind 
the intersection point due to the continuous distribution of sources 
within the interval may be obtained from equation (l) by use of appropriate 
limits of integration. The variation of pressure with x given by this 
expression is closely approximated by equation (A2) for the source 
location just determined. v - 

It is therefore evident that the pressure on the surface of a cone 
at and behind the intersection point due to a continuous distribution 
of sources along the axis within an interval Ax from the vertex is 
closely approximated by a single source located Ax/4 ahead of the end 
of the interval. Further computation along these lines has indicated 
that this principle may be applied to an interval located in any part 
of a conical or parabolic body of revolution with reasonable accuracy, 
provided that the pressure function Q is computed (equations (2) and (4)) 
from the radius and slope of the body surface at the source location. 

It appears reasonable to assume, therefore, that within the general 
limitations of the linear theory, the principle may be applied to 
arbitrary bodies. 
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TABLE I 

YALUES FOB GRID (FIG- l) OF LIKES OF CONSTANT ERASURE COEFFICIEIIT 
IK x,Pr ELANE CALCTILAIEED FROM EQUATION (3); SOURCE AT ORIGIN; 
ERESSUKE FUNCTION Q = O.Ol 


\pr 

p \ 

0 

\ 

0.01 

0.02 

0.03 

0.05 

0.09 


X 

03 

0 

0.0100 

0.0200 

0.0300 

0.0500 

0.0900 

2.00 

.0100 

.0141 

.0224 

.0316 

.0510 

.0908 

1.00 

.0200 

.0224 

.0283 

.0361 

.0538 

.0922 

.70 

.0286 

.0303 

.0349 

.0414 

.0576 

.0944 

.50 

.0400 

.0412 

.0447 

.0500 

.o 64 o 

.0995 

• 30 

.0667 

.0674 

.0696 

.0731 

.0834 

.1120 

.20 

.1000 

.1005 

.1020 

.1044 

.1118 

.1345 

.16 

.1250 

.1254 

.1266 

.1286 

. 13^6 

.1540 

.10 

.2000 

.2002 

.2010 

.2022 

.2062 

.2193 

.08 

.2500 

.2502 

.2508 

.2518 

.2550 

.2657 

.05 

.4000 

.4001 

.4005 

.4011 

.4031 

.4100 

• 03 

.6667 

.6667 

.6670 

.6674 

.6685 

.6727 

.02 

1.0000 

1.0000 

1.0002 

1.0004 

1.0012 

1.0040 
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T A~ffr.lt XT 

FORM FOB TABOLAimt} TALOES OF B OH ST3HFACE OS' BOOT EtE TO EACH SOURCE AHD -VALTIES OF I AT MICH 
TEET OCCUR. EXAHPTE: EARABGLIC BOOT OF FTHKWSB RATIO 12 AT A MACH HDMHER OF 1.4l 

[Ondorlined. and toacEated values are far erample] . ' 



-0.4625 

—4125 

-.3625 

-•3125 

[-.2625] 

-.2125 
-.1625 
-.n 25 
-.0625 
—0125 
.0375 
.0875 

.1375 

.1875 

.2375 

.2875 

•3375 

.3875 

•4375 

.4875 


O.IO85 
.0725 
.0398 
.0112 
» -.0118 

» -.0315 

“ -.0475 

» —0594 

“ —.0660 

» -.0692 

CO —.0683 

co -.0633 

CO -.0537 

co -.0403 

CO -.0222 

CO 0 

CO .0255 

CO .0553 

05 .0900 

CO .1290 


0.0542 
.0362 
.0199 
.0056 
£-.00 
-.01 
-.0238 

-.0297 

-.0330 

-.0346 

—0342 

-.0316 

-.0268 

-.0202 

-.0111 

0 

.0128 

-0276 

.0450 


0.0217 
.0145 
.0080 
.0022 
-.0024 
— 0063 
-.0095 
-.0119 
-.0132 
-.0138 
-0137 
-.0127 
-.0107 
-.0081 
-.0044 
0 

.0051 

.0111 


0.0174 0.01 
.0116 
.0064 
.0018 
[-0019] 

-.0050 
—0076 
—0095 
-.0106 
—0111 -.00 
-.0109 -.00 
—aioi —0063 
—0086 —0054 
— 0064 — 0040 
—0036 
0 

.0041 


0.0054 0.0033 
.0036 .0022 
.0020 .0012 
.0006 .0003 
—0006 —0004 
-.0016 -.0009 
—0024 
—0030 
-0033 
-.0035 
—0034 
—.0032 



0.1085 

.0725 

..*0398 

.0112 

£—0118] 

• -.0315 
—0475 
- -10594 
—0660 
—0692 
-.0683 
-.0633 
-•0537 

-.0403 

■-.0222 

0 

.0255 

•0553 

..0900 

.1290 


n,NACA^- 







































X(Qrtd) 

; . , l 

' ’ ' ■' , . V . 

Figure 2.- Parabolic "body of fineness ratio 12 plotted on x,Pr plane for M = l,4l. Grid (fig. 1) 

is shown superimposed .on "body with origin of grid located at source location £ = -0.2625. 
Ordinate scale is 10 times abscissa scale . 
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Figure 4.- Variation of pressure coefficient P along the surface of 
the parabolic body of fineness ratio 12 due to each source (see 
table H) at a Mach number of lAl. 






Figure 5 .- Comparison of pressure distributions 

-T» — _*U - T 1 - 1.. J— __C» J31— - — .XJ _ irt 

XUX' tl LKXUJ UX X XXLtJXLtJbb X'tL U1U J-C. 

computed at a Mach number of lXl Py tile 
point -source method and Py the method of 
Jones and Margo He (reference 4) . 



-S -A r3 -Z V 0 J 2 2 A .5 

%■ 

Figure 6.- Comparison of pressure distriPutions 
for a parabolic PodLy of fineness ratio 6 
computed at a Mach number of lXl Py the 
point -source method and Py the method of 
Jones and Margo lis (reference 4) . 
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Figure 7.- Comparison of calculated, and experi- 
mental pressure distributions for the 
A tV IP at a Mach number of I.87. Experi- 
mental points and solid fairing tafcen 
:from reference 7. 
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Figure 8.- Comparison of pressure distributions 
calculated by the present method and by the 
characteristics methods of Sauer and Tollmein 
with vind -tunnel results for the AtVIP at a 
Mach number of 1.87. 


KA.CA TS Ho. 1768 




ro 

■p- 





Figure 9 .- Comparison of pressure distributions Figure 10.- Comparison of pressure distributions 

calculated by the present method and by the calculated at a Mach number of 1.41 by the 

characteristics methods of Sauer and Tollmein point -source method for a parabolic body of 

vith vind -tunnel results for the A4Y1P at a fineness ratio 12 and for the free-fall test 

Mach number of 1 . 56 . body of fineness ratio 12. 
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